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For Mars, 57,633 craters from the manually assembled catalogues and 72,668 additional craters identified 
using several crater detection algorithms (CDAs) have been merged into the MAI 30301 GT catalogue. 
By contrast, for the Moon the most complete previous catalogue contains only 14,923 craters. Two recent 
missions provided higher-quality digital elevation maps (DEMs): SELENE (in 1/16° resolution) and Lunar 
Reconnaissance Orbiter (we used up to 1/512°). This was the main motivation for work on the new Crater 
Shape-based interpolation module, which improves previous CDA as follows: (1 ) it decreases the number of 
false-detections for the required number of true detections; (2) it improves detection capabilities for very 
small craters; and (3) it provides more accurate automated measurements of craters’ properties. The results 
are: (1) LU60645GT, which is currently the most complete (up to ~D > 8 km) catalogue of Lunar craters; 
and (2) MA132843GT catalogue of Martian craters complete up to > 2 km, which is the extension of the 
previous MAI 30301 GT catalogue. As previously achieved for Mars, LU60645GT provides all properties that 
were provided by the previous Lunar catalogues, plus: (1) correlation between morphological descriptors 
from used catalogues; (2) correlation between manually assigned attributes and automated measurements; 
(3) average errors and their standard deviations for manually and automatically assigned attributes such as 
position coordinates, diameter, depth/diameter ratio, etc; and (4) a review of positional accuracy of used 
datasets. Additionally, surface dating could potentially be improved with the exhaustiveness of this new 
catalogue. The accompanying results are: (1) the possibility of comparing a large number of Lunar and 
Martian craters, of e.g. depth/diameter ratio and 2D profiles; (2) utilisation of a method for re-projection of 
datasets and catalogues, which is very useful for craters that are very close to poles; and (3) the extension of 
the previous framework for evaluation of CDAs with datasets and ground-truth catalogue for the Moon. 

© 201 1 Elsevier Ltd. All rights reserved. 


1. Introduction 

Recently, all the craters from the major currently available 
manually assembled catalogues have been merged into the 
MA57633GT catalogue with 57,633 known Martian impact-craters 
(Salamuniccar and Loncaric, 2008b). In addition, using our crater 
detection algorithm (CDA) and 1/128° MOLA data, 57,592 previously 
uncatalogued craters have been identified resulting in the 
MA115225GT catalogue (Salamuniccar and Loncaric, 2010a). This 
catalogue has been additionally extended using: (1) MA75919T 
catalogue, which is the result of the CDA from Stepinski and 
Urbach (2008); (2) the CDA from our previous work (Salamuniccar 
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and Loncaric, 2010a), wherein instead of the Canny we used Shen- 
Castan edge detector; and (3) the CDA for optical images (Bandeira 
et al., 2007) for four selected regions of interest. This resulted in 
MAI 30301 GT catalogue (Salamuniccar et al., 2011b), which is at the 
time of writing the most complete previously available public 
catalogue of Martian craters. By contrast, for the Moon the most 
complete previous catalogue contains only 14,923 craters 
(Rodionova et al., 1987). The size difference between Mars and the 
Moon only partially explains the difference, because on the Moon the 
average density of the largest craters is greater than on Mars 
according to size-frequency distribution of the largest craters, for 
which these catalogues are mostly complete (see below; Section 3.2 
- the distribution of craters from corresponding catalogues). Two 
recent missions provided higher-quality digital elevation maps 
(DEMs): (1) SELENE laser altimetry (LALT) dataset in 1/16° resolution 
(Araki et al., 2009); and (2) Lunar Reconnaissance Orbiter ( LRO ) Lunar 
orbiter laser altimeter (LOLA) dataset (Smith et al., 2010), wherein 
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we used resolutions up to 1 / 512°. This was a motivation for our work 
on a considerably more complete catalogue of Lunar craters. 

Research on CDAs is challenging for numerous reasons. On 
Mars, there are many different crater shapes depending on their 
interior morphologies (central peaks, peak rings, central pits, and 
wall terraces) and ejecta structures (pedestal, pancake, rampart, 
lobate, fluidized, radial or lunar-like, transitional or diverse) 
(Barlow et al., 2003). There is also a simple-complex transition 
from the smaller, mostly very circular bowl-shaped craters, to 
larger complex craters with central peaks, and to the largest 
multi-ring impact basins (Melosh and Ivanov, 1999). In general, 
CDAs are based on a large number of methods, including circle/ 
ellipse detection (Cooper, 2003; Cooper and Cowan, 2004; Flores- 
Mendez and Suarez-Cervantes, 2009; Krogli and Dypvik, 2010), 
probability volume created by template matching (Bandeira et al., 
2007), machine-learning (Stepinski and Urbach, 2008), etc. 
An overview of 112 publications related to CDAs has been 
published in two recent papers (Salamuniccar and Loncaric, 
2008a; Salamuniccar et al., 2011b). The list is constantly increas- 
ing, as is e.g. the work regarding feature selection and boosting 
(Ding et al., 2010) and entropic quadtrees (Vetro and Simovici, 
2010). However, the state of the art in image-analysis/object- 
recognition still does not offer an answer on how to create a CDA 
that is as robust as the scientific community would desire. Hence, 
it is not surprising that several research groups are working in 
this still young field, trying to achieve this goal. Several CDAs can 
process global datasets (Michael, 2003; Salamuniccar and 
Loncaric, 2010a; Stepinski and Urbach, 2008). The experience 
with our CDA is as follows: (1 ) it is necessary to manually process 
a large number of crater-candidates in order to considerably 
extend ground-truth (GT) catalogue; (2) improvements of detec- 
tion capabilities for very small craters increase the number of 
detected craters from any given dataset; and (3) the inaccuracy in 
measurement of craters’ properties including 2D profiles increases 
with the decrease in their size. This was a motivation for the work 
on the new Crater Shape-based interpolation, which considerably 
improves the previous CDA in addressing these issues. 

The above led to the concurrent work on new catalogues and 
improvements of the previous CDA (Salamuniccar and Loncaric, 
2010a). Resulting catalogues are: (1) MA132843GT for Mars; and 
(2) LU60645GT for the Moon, created using LALT (Salamuniccar and 
Loncaric, 2010b) and LOLA (Salamuniccar et al., 2011a) datasets. 

The rest of this paper is organised as follows. In Section 2, the 
methods and datasets of this paper are presented in detail, as well as 
used and assembled datasets. The results are presented in Section 3, 
and the conclusion is given in Section 4. Readers mostly interested 
into the improved CDA from this paper should be primarily focused 
to Sections 2.1, 2.2, and 3.1, readers mostly interested into the 
resulting catalogues should be primarily focused to Sections 2.3, 2.4, 
and 3.2., while Section 3.3 demonstrates how new crater catalogues 
can be used in combination with the improved CDA which provides 
automated measurements, in order to outline differences between 
Martian and Lunar craters. 


2. Methods and datasets 

The developed methods and the datasets used are described in 
the following order: 

• improvement of crater detection algorithm using new Crater 
Shape-based interpolation; 

• utilisation of method for re-projection of datasets and 
catalogues; 

• datasets used for crater detection and evaluation; and 

• previous catalogues of Lunar craters. 


The following are made available online (http://informatika. 
tvz.hr/index.php?pred=17461): the source code (Craters5_77.zip) 
of the interpolation-based CDA presented in this paper; the 
Topolyzer application, which implements methods for the eva- 
luation of CDAs and the registration to ground-truth (GT) catalo- 
gues; and the catalogues LU60645GT and MA132843GT (as well as 
previous versions). In addition, the various datasets are available on 
request, in order to provide reproducibility of the methods pre- 
sented in this paper and inter-comparison with other CDAs. 

2.1. Improvement of crater detection algorithm using new Crater 
Shape-based interpolation 

The CDA from this paper, which uses multi-resolution image 
analysis and the Crater Shape-based interpolation algorithm 
described in following paragraphs, is shown in Fig. 1. The algorithm 
is an extension of our previous work (Salamuniccar and Loncaric, 



Fig. 1 . Schematic flow diagram of overall organisation and processing steps for 
our CDA. 
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2010a). The algorithm utilises fuzzy edge detection to extract crater 
features required for the fuzzy Radon-Hough (RH) transform search 
for maxima in the parameter space. Following this step, DEM and 
parameter-space data are interpolated at run-time in order to be 
used for morphometry measurements (of depth/diameter ratio, 
circularity, topographic-cross-profile, rim, central peak, and radial 
range where the crater is preserved) and parameter space analysis 
(the circularity of votes is higher at craters’ centres than in the 
centres of false detections). The interpolation at run-time is used to 
dynamically calculate interpolated values on demand, therefore 
reducing memory requirements. Next, the previously computed 
probability that a detected feature is a crater is modified based on 
its associated diameter-range (not compared to the resolution of the 
image currently used in the multi-resolution setup, but according to 
classification into one of diameter-ranges - see below; Section 3.2 - 
the distribution of craters from corresponding catalogues). These 
steps are repeated for all radii and for all data resolutions. Finally, 
we perform slip-tuning of the crater-candidates, and removal of 
multiple detections. Because the previous version of CDA 
(Salamuniccar and Loncaric, 2010a) was used only for Mars, the 
planetary-radius had been hard-coded. In the new version, it 
becomes one parameter, so that the same CDA can be used for 
Mars, the Moon, and potentially for other planetary bodies. This 


value is used inside the CDA for the computation of slopes and 
diameter-range of crater according to its size. The proposed CDA can 
use various edge detection methods such as Canny (1986) and Shen 
and Castan (1992). The experiments conducted in this work are 
based on the Canny (1986) edge detection method, and have shown 
considerably improved crater detection. 

The proposed DEM-based CDA is based on our previous CDA 
(Salamuniccar and Loncaric, 2010a). It has been improved using a 
specially developed method for interpolation of DEM and para- 
meter-space data, which is suitable for detection of small craters. 
The most commonly-used interpolating functions are: nearest 
neighbour, bilinear, bicubic, quadratic splines, cubic B-splines, 
higher-order B-splines, Catmull-Rom cardinal splines, Gaussians, 
and truncated sine (Zitova and Flusser, 2003). However, in the case 
of their application to CD As, we should consider that: (1) the 
optimal implementation of interpolation should be computed in 
runtime only for required DEM coordinates, as opposed to a static 
interpolation of the whole data set (this reduces computational 
complexity and memory requirements); (2) the interpolation itself 
should not introduce artifacts that will decrease the accuracy or 
precision of the measurements obtained via CDA (e.g. compiled 
depth/diameter values); and (3) for the smallest detectable craters 
we have sparse sampling, wherein each crater is represented with 
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Fig. 2. Illustration of the Crater Shape-based interpolation used for the CDA from this paper: (a) the 2D topographic profile of crater; (b) the representation with nearest- 
neighbour interpolation of the elevation; (c) the simple linear interpolation; (d) the second-order interpolation based on the additional interpolation (see the bottom 
frame), but only for distances smaller than 1 pixel; and (e) the comparisons which show that (d) is a better approximation of the profile shown in (a) than the initial case 
(b) and previous case (c). 
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a small number of pixels. On the other hand, the importance of good 
interpolation lies in: (1) better detection capabilities, particularly for 
small craters; and (2) more accurate measurements of crater 
properties. Based on these considerations, we propose the following 
solution. 

Crater Shape-based interpolation used for our CDA is illu- 
strated in Fig. 2. Frame (a) shows the topographic profile of a 
crater from MA132843GT. Frame (b) shows a pixel representation 
(known also as the nearest neighbour function), whose pixel size 
with artificially degraded resolution illustrates the effect of the 
digitalization of input values. It also enables us to see how various 
interpolation methods compare to the original full-resolution 
data. The result of simple linear interpolation is shown in frame 
(c). It shows an improvement for most segments, except very 
close to the crater centre. In such cases there is a decrease of 
accuracy (the curve shown in frame (c) has a first derivative 
discontinuity at the crater centre, which is not the case for the 
curve shown in frame (a)). In order to overcome this interpolation 
artifact, we modify the interpolation between the crater centre 
and the nearest pixels, as shown in frame (d). The additional 
interpolation is shown in the bottom frame of Fig. 2. It assumes a 
parabolic crater shape, which is the case for the majority of small 
(bowl-shaped) craters, and which is also appropriate for larger 
craters with central peaks or pits (concave vs. convex parabola). 
As shown in frame (e), this approximates original (full-resolution) 
values in all segments considerably better than previous cases 
(b) and (c). At the same time, approach (d) is computationally 
very efficient, because simple linear interpolations are computed 
in only two stages, the second of which is only executed within a 
one-pixel distance to the crater centre. Even within a one-pixel 
distance to the crater centre, this is considerably faster than 
alternative, more complex interpolation approaches (e.g. the 


computation of coefficients of higher degree polynomials). There- 
fore, we selected approach (d) for our CDA. 

2.2. Utilisation of method for re-projection of datasets and 
catalogues 

In previous work (Salamuniccar et al., 2011b), we provided an 
orthographic projection view with the crater of interest in the 
centre for better evaluation capabilities of northern-most and 
southern-most craters. This approach was sufficient for Mars 
where all the craters from MAI 30301 GT (Salamuniccar et al., 
2011b) are located between 87°N and 87°S. Flowever, this is not 
the case for the Moon, where there are no polar-layer-deposits, 
and where craters exist at the most extreme polar latitudes. In 
addition to orthographic projections, we perform a rotation of the 
entire DEM and optical global datasets, such that the geographical 
North and South Poles are rotated on the equator of the new 
coordinate system. The method works as follows: (1) a global 
dataset is re-projected from the simple cylindrical projection to 
the sphere; (2) the sphere is rotated by a parameterized amount 
around the x, y, and z axes; (3) the global dataset is re-projected 
back from the sphere to the simple cylindrical projection. 
The same method can be used for crater locations, in both direc- 
tions: (1) to convert geographical crater coordinates to comply with 
the re-projected datasets; and (2) to convert the rotated coordinates 
back into the original form, to comply with the initial datasets. 

2.3. Datasets used for crater detection and evaluation 

For Mars, we used datasets prepared in our previous work 
(Salamuniccar et al., 2011b). For the Moon, we used the 1/16° 
SELENE LALT (Araki et al., 2009) and up to 1/512° LRO LOLA (Smith 



Fig. 3. As shown for one smaller (D = 5.1 km) V-shaped crater (A-l ) and one slightly larger (D = 7.8 km) bowl-shaped crater with an approximately parabolic interior profile 
(B-l), the increase of the number of radial-resolution scans cannot solve the pixelization-side-effect problem (A-3 compared to A-2, B-3 compared to B-2). In contrast, the 
Crater Shape-based interpolation from this paper successfully solves this problem (A-4 compared to A-3 and A-2, B-4 compared to B-3 and B-2). The resulting crater 3D 
shape and 2D topographic profiles (A-4 and B-4) show no artifact resulting from the sub-optimal data-processing. 
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et al, 2010) global DEMs, digitised at 1-m vertical resolution. 
Unlike previous laser altimeters, which mapped planetary bodies, 
such as Mars orbiter laser altimeter ( MOLA ), LOLA is a five-beam 
laser altimeter with a firing rate of 28 Hz, providing up to 140 
measurements per second. The five-beam pattern was designed 
to characterize the lunar surface at scales relevant to human and 
robotic exploration, and it translates to an unprecedented along- 
track measurement spacing of ~10m (compared to about 1.7 km 
for SELENE). LOLA continues to operate, improving its coverage 
(especially, reducing its average longitudinal gap between tracks). 
The global Lunar optical image mosaics used during the registra- 
tion of new craters in the GT catalogue are: 1/64° LASGW (Lunar 
Airbrushed Shaded Relief Warped to ULCN2005); 1/256° CLEM- 
BASE (Clementine Basemap Mosaic version 2); and 1/512° LOM 
(Lunar Orbiter Mosaic). In the case of LALT and LOLA DEM datasets 
and LASGW optical dataset, post-processing tasks includes only 
down-sampling and up-sampling to other resolutions, and ren- 
dering of optical BMP files. For CLEMBASE, we performed an 
interpolation fix (as done for Martian THEMIS-DIR ) and bright- 
ness/contrast correction fixes (as done for Martian MOC). The LOM 
dataset was optimised for memory usage purposes. 

2.4. Previous catalogues of Lunar craters 

Two Lunar catalogues developed by previous researchers are: 
(1) McDowell (2007) catalogue, which contains 8639 (named) 
craters (hereafter LU8639N); and (2) Rodionova et al. (1987) 
catalogue (hereafter LU14923R), which is the most complete 
catalogue of Lunar craters prior to this work, containing 14,923 
craters (1394 identified by name). An additional Lunar catalogue 
(Head et al., 2010; Radish et al., 2011), which contains 5185 
craters (hereafter LU5185H) was developed very recently, and 
was not available when we had completed LU58357GT released 
earlier in order to facilitate the use of LRO data to survey 
and estimate impact melt volumes in small- to medium-sized 
craters (Mazarico et al., 2011). LU5185H was thus included in 
LU60645GT. 


3. Results and discussion 

In order to evaluate the performance of the proposed CDA, 
according to the framework for evaluation of CD As (Salamuniccar 
et al., 2011b), the resulting new catalogues MA132843GT and 
LU60645GT are used in addition to the previously assembled 
catalogue MA57633GT (Salamuniccar and Loncaric, 2008b). 
The results and discussion are given in the following order: 

• improved crater detection algorithm based on new Crater 
Shape-based interpolation; 

• assembled new catalogues MA132843GT and LU60645GT of 
Martian and Lunar craters; and 

• differences between Martian and Lunar craters in depth/ 
diameter ratio and 2D profiles. 

Section 3.1 additionally includes: (1) interpolation-based 
improvements in 3D-shapes of craters and 2D-topography-pro- 
files; and (2) evaluation of achieved performance of interpolation- 
based crater detection algorithm. Section 3.2 additionally 
includes: (1) verification of assigned coordinates and diameters 
to craters from polar- regions; (2) typical newly catalogued craters 
from MA132843GT and LU60645GT catalogues; (3) statistics of 
LU60645GT and MA132843GT catalogues; and (4) advances in 
integration of diverse morphological descriptors from previous 
catalogues. Section 3.3 describes one possible application of crater 


catalogues - insight into differences between Martian and Lunar 
craters. 


3 A. Improved crater detection algorithm based on new Crater 
Shape-based interpolation 

The Crater Shape-based interpolation used within the CDA 
results in improved crater 3D-shapes and 2D topographic profiles, 
as shown in Fig. 3 for two fresh craters with large depth/diameter 
ratios. The curve shown in Fig. 2 frame (c) has a discontinuation of 
first derivation at the crater rim as well (and not only at the crater 
centre as previously discussed), as well as the curve shown in frame 
(d), while for the curve shown in frame (a) this is not the case. From 
this perspective, another important observation is that the disconti- 
nuity of the first derivative at the crater rim does not appear in 
Fig. 3. This is so because this problem has been solved by averaging 



Fig. 4. Part of the 1/128° LOLA dataset and output of our Crater Shape-based 
interpolation CDA within (blue) and outside (red, not shown in the bottom frame) 
the selected region (top); crater-candidates with the assigned probability larger than 
the used threshold which defines those manually evaluated, true (blue) and false 
(red, labelled with white arrow-pointers) detections, according to the automated 
verification using LU60645GT (bottom). (For interpretation of the references to color 
in this figure legend, the reader is referred to the web version of this article). 


Table 1 

Definitions and graphs (TP=True Positives; FP= False Positives; FN = False Nega- 
tives; GT= Ground Truth; TDR=True Detection Rate; FDR= False Detection Rate; 
ROC=Receiver Operating Characteristics; F-ROC=Free-response ROC; ROC’ = the 
closest approximation that can be achieved for ROC; Q= Quality percentage) used 
for evaluation of CDAs (Salamuniccar et al., 2011b). 


Used definitions 

Used graphs 


(horizontal range/vertical range) 

GT=TP + FN 

F-ROC [TP/FP] 

TDR=TP/(TP + FN)=TP/GT 

ROC’ [TDR/FDR] 

FDR=FP/(TP + FP) 

Q4Q/FDR] 

Q=TP/(TP+FP + FN)=TP/( GT + FP ) 
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a large number of profiles (we average 176 profiles at different 
azimuths/angles of the same crater, partially because our CDA 
searches for craters in radius range between 2 and 28 pixels and 
in such case the length of rim in pixels is 28 x 2n= ~176). 

The difference compared to the CDA from previous work 
(Salamuniccar and Loncaric, 2010a), which utilised only the Martian 
DEM ( MOLA ) is that in this work the Lunar DEMs have been used 
(LALT, LOLA ) as well. The proposed interpolation-based CDA can 
successfully detect Lunar craters as shown in Fig. 4: (1) most of the 
craters are successfully detected; (2) the number of false-detections is 
low; and (3) coordinates and diameters are precisely defined 
for a large majority of the craters. As shown, for this region 
(xl = “ - 96.578 132” yl = “20.1 953 15” x2 = “- 90.929678” y2 = 
“16.921865”; xl and x2 are longitudes, yl and y2 are latitudes of 
the corners, in the east/planetocentric coordinate system) there are 
17 true detections and 6 false detections. All false detections have a 
diameter smaller than 3.317 km and for larger craters there are no 
false detections within this region. For parameterization of the 
interpolation-based CDA proposed with this paper, we used the 
planetary radius value for Mars of 3390 km, and for the Moon 
1737.4 km. 

The evaluation was performed using the framework for eva- 
luation of CDAs (for definition see Table 1) from previous work 


(Salamuniccar et al., 2011b), using F-ROC (Salamuniccar and 
Loncaric, 2008a) and ROC’ (Salamuniccar et al., 2011b). The ROC’ 
evaluations for Mars show that the proposed CDA is better than 
all other available CDAs, according to the area under the ROC’ 
curve criterion (AUROC’ is 50.1% for the CDA by Stepinski and 
Urbach (2008), 54.0% for the CDA by Salamuniccar and Loncaric 
(2010a), and 57.4% for the proposed interpolation-based CDA). 
The ROC’ evaluations based on lunar datasets (shown in Fig. 5 
bottom) show that the proposed interpolation-based CDA is also 
better for Lunar datasets than the CDA from the previous work 
(Salamuniccar and Loncaric, 2010a), while other CDAs have not 
been used before for global processing of Lunar DEM. 

In order to compare the relative improvements of overall 
achieved performance between our previous CDA (Salamuniccar 
and Loncaric, 2010a) and the CDA from this paper (wherein the 
only difference is in the usage of Crater Shape-based interpola- 
tion) with the improvements between the CDA from Stepinski 
and Urbach (2008) and our previous CDA (Salamuniccar and 
Loncaric, 2010a), the following method was used. Let us suppose 
that CDA-A, CDA-B, and CDA-C achieved performance A, B and C, 
wherein: (1 ) each value is in the normalised range between 0 and 
100%; (2) larger value means better performance; and (3) A<B 
and B<C. In this case, for CDA-B possible improvement is (1 - A) 
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Fig. 5. According to the F-ROC evaluations for Mars (top-left) CDA #2 is superior to CDA #1 (Stepinski and Urbach, 2008), while CDA #3 provides better results than the 
other two. According to the ROC evaluations for Mars (top-right) CDA #3 is superior to CDA #2, but CDA #1 shows better results for false detection rates between 5% and 
30%. According to the ROC’ evaluations for the Moon (bottom-left and bottom-right), CDA #3 is superior to CDA #2 like in the case for Mars. For 1/128° LOLA data obtained 
results are better than for 1/16° LALT or 1/64° LOLA data. 
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and achieved improvement is ( B-A ), while for CDA-C possible 
improvement is (1 -B) while achieved improvement is (C-B). 
Using this, the ratio between achieved and possible improvement 
is R b =(B-A)/( 1 -A) for CDA-B and R c =(C-B)/( 1 -B) for CDA-C. 
At the end, these two values are compared, having J B =R B - 100/ 
(R b -\-Rc) and I c =Rc- 100/(R B +R C )- For AUROC values from the 
previous paragraph, computed values are 50.6% and 49.4%. 
Similarly, we can compute for AUFROCi x , in which case we are 
not interested in areas under F-ROC for a horizontal range larger 
than It (t is the number of craters in GT catalogue) because in 
such a case the number of false detections is always larger than 
the number of correct detections. In order to take into account the 
fact that some craters from MA132843GT catalogue cannot be 


detected by either CDA because of dataset limitations, an addi- 
tional GT catalogue was created from which 19,886 craters not 
detectable by either CDA have been removed. The overall results 
are shown in Fig. 6, wherein we also averaged the results from 
individual computations in order to have as realistic as possible a 
comparison of achievable performances between different CDAs. 

3.2. Assembled new catalogues MA132843GT and LU60645GT of 
Martian and Lunar craters 

The results of the method for re-projection of datasets and 
catalogues are shown in Fig. 7. Using this method, the coordinates 
and diameters of 22 Lunar craters were corrected. After that, 
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Fig. 6. Comparison of the relative improvements of overall achieved performance between: (1) the CDA from Stepinski and Urbach (2008) and our previous CDA 
(Salamuniccar and Loncaric, 2010a) - left parts (blue); and (2) our previous CDA (Salamuniccar and Loncaric, 2010a) and the CDA from this paper wherein the only 
difference, which influences performance is in the usage of Crater Shape-based interpolation - right parts (green). (For interpretation of the references to color in this figure 
legend, the reader is referred to the web version of this article). 
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Fig. 7. Standard (A-1 and B-1) and rotated (A-2 and B-2) view for 1/1° MOLA (A-1 and A-2) and 1/1° LOLA (B-1 and B-2) data. The rotated view was used to check the 
accuracy of assigned coordinates and diameters of craters, as shown for 1/64° LOLA (C-l and C-2) North-Pole (C-l) and South-Pole (C-2). 
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we checked the 50 northern-most and 50 southern-most craters, 
in order to confirm that additional corrections are not required for 
LU60645GT. A similar check was performed for MA132843GT, 
wherein corrections have not been required. 

To extend the GT catalogues, we used the framework for 
evaluation of CD As from previous work (Salamuniccar et al., 
2011b), where all crater-candidates proposed by our CDA are 
manually evaluated for rejection, or correction of coordinates and 
diameter and inclusion into the GT catalogues. The main phases of 
the new work on GT catalogues for Mars and the Moon are shown 
in Fig. 8. We rejected some of the crater-candidates before 
manual evaluation, using the framework for evaluation of CDAs 
from previous work (Salamuniccar et al., 2011b). In the first step 
(Mars - 1/128° MOLA), we found the maximum of Q (Quality 
percentage) curve (for definition see Table 1) using MAI 30301 GT 
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Fig. 8. Main phases of the new work on GT catalogues for Mars (top) and the 
Moon (bottom) in which the new work is based on the interpolation-based CDA. 


and the associated probability threshold ~ 0.3907 (hereafter t) 
and removed all crater-candidates with associated probabilities 
smaller than t. The recomputed t for MA132843GT is the same 
number. The value of t computed for Mars has been used for 
Lunar catalogues as well because it was obtained using GT 
catalogue, which was much more complete than that available 
for the Moon at that time. When we assembled LU60645GT we 
recomputed the Lunar value of t. It is remarkably similar to 
the previously computed for Mars (difference is only of the order 
10 -5 ). This is explainable by completeness of LU60645GT in 
comparison with MA132843GT and the robustness of the inter- 
polation-based CDA proposed by this paper, regarding detection 
of craters on two different planetary bodies (Moon and Mars). 

As shown in Fig. 9 (top), typical new entries in Martian GT 
catalogue are mostly small craters, and some larger much 
degraded craters. These are the craters that the new interpola- 
tion-based CDA from this paper successfully detected, which was 
not the case for our previous CDA (Salamuniccar and Loncaric, 
2010a). As shown in Fig. 9 (bottom), there are many more new 
entries in Lunar GT catalogue than in the Martian one. The reason 
is that preceding Lunar catalogues were considerably less com- 
plete than the preceding MAI 30301 GT Martian catalogue. 

The distribution of craters from corresponding catalogues and 
classification into one of diameter-ranges is given in Fig. 10. The 
completeness is estimated using the Craterstats (Michael, 2010) 
programme. The flattening of lines at small diameters is related to 
the input dataset resolution. The CDA detected numerous small 
craters (D < 8 km) in the LOLA dataset, which could not be verified 
due to inadequate high-resolution imagery coverage, and there- 
fore has not been included into LU60645GT. In addition, we 
expect a correlation between poorer imagery coverage and 
reduced density of small craters (D < 8 km) (e.g. fewer smaller 
detected craters validated on the farside or higher-latitudes, 
where imagery coverage was worse). 

The catalogue integration is useful for researchers interested in 
craters on Mars and the Moon. The new catalogue LU60645GT is a 
superset of the previous Lunar catalogues, as shown in Fig. 11. 
With our approach, LU60645GT contains everything that was 
included in these catalogues, plus: (1) the correlation between 
various morphological descriptors from the catalogues used; 
(2) the correlation between manually assigned attributes and 
automated depth/diameter measurements from LU60645GT by 
our CDA; (3) average errors and their standard deviations for 
manually and automatically assigned attributes, such as position 
coordinates, diameter, depth/diameter ratio, etc; and (4) posi- 
tional accuracy of features in the used datasets, wherein our 
catalogue contains 60,645 cross-references between each of 
the used datasets. Additionally, surface dating could potentially 
be improved with the exhaustiveness of this new catalogue. 
The same kind of contribution was also made in the case of the 
new Martian catalogue MA132843GT, as described into more detail 
in the previous work for MAI 30301 GT (Salamuniccar et al., 2011b). 

3.3. Differences between Martian and Lunar craters in depth/ 
diameter ratio and 2D profiles 

The Crater Shape-based interpolation CDA from this paper can be 
used in combination with the new LU60645GT and MA132843GT 
catalogues of Lunar and Martian impact craters in order to provide 
automated morphometry measurements of depth/diameter ratio, 
topographic-cross-profiles, etc. (Salamuniccar and Loncaric, 2009, 
2010a). Now that catalogue LU60645GT exists, in addition to 
MA132843GT (extension of the previously available MA130301GT), 
it is possible to compare morphometry measurements of Martian 
and Lunar impact craters. As shown in Fig. 12 (top), the transition 
of the depth/diameter ratio of the youngest craters on Mars occurs 




244 


G. Salamuniccar et al. / Planetary and Space Science 60 (2012) 236-247 



Fig. 9. Previous MA130301GT catalogue (green) and additions (yellow, labelled with white arrow-pointers) of the new MA132843GT catalogue of Martian craters (top). 
The catalogue with 18,229 craters assembled by human labelers (green) and additions (yellow, labelled with white arrow- pointers) of the new LU60645GT catalogue of 
Lunar craters (bottom). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article). 
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Fig. 11. The resulting LU60645GT catalogue contains LU8639N, LU14923R, and LU5185H, wherein it is possible to compare the properties of numerous craters between 
catalogues. 


at D=5.8 km, while on the Moon it occurs at D= 18.4 km. According 
to its diameter, we can classify each crater into one of the seven 
following ranges: (1) lkm<D<2km; (2) 2km<D<4km; 
(3) 4 km < D < 8 km; (4) 8 km < D < 1 6 km; (5)16 km < D < 32 km; 
(6) 32km<D<64km; and (7) 64km<D. As shown in Fig. 12 
(centre), it is possible to compare average topography-cross-profiles 
of Lunar and Martian craters. They follow the trend of the youngest 
craters, wherein for Lunar craters depth/diameter ratio and height of 
crater rim is approximately two times larger than for Martian 
craters. For each average topographic-cross-profile, using d/D ratio, 
all craters from a selected range are in addition classified according 
to whether they have higher or lower d/D ratio. The whole 
procedure is recursively repeated 3 times. The result is that the 
craters within each diameter range are additionally classified into 
8 sub-groups according to their d/D value. As shown in Fig. 12 
(bottom), the above trend is consistent for all ranges from the 
freshest to the most heavily degraded/eroded craters. 


4. Conclusions 

The main results of this paper are: (1) the new catalogue 
LU60645GT of Lunar craters (complete up to~D>8km), 
which contains a considerably larger number of craters and 
associated attributes than any other currently available catalogue 
of Lunar craters; (2) the new catalogue MA132843GT of Martian 
craters (an expansion of the previous MAI 30301 GT), which is the 
most complete (up to~D>2 km) publicly available catalogue of 
Martian craters; (3) considerable improvement of previous CDA 
(Salamuniccar and Loncaric, 2010a) using new Crater Shape- 
based interpolation method from this paper, wherein the 
improved CDA was successfully used during work on LU60645GT 
and MA132843GT. As previously achieved for Mars with our 
previous catalogue MAI 30301 GT (Salamuniccar et al., 2011b) 
and catalogue MA132843GT from this paper, for the Moon the 


new LU60645GT catalogue provides everything that was provided 
by the previous Lunar catalogues, plus numerous new possibili- 
ties as enumerated in Section 3.2. The accompanying results of 
this paper are: (1) the possibility of comparing Lunar and Martian 
craters regarding depth/diameter ratio, 2D profiles, etc., wherein a 
considerably larger number of craters than available from the 
previous work leads to globally more representative, statistically 
more significant, more precise results; (2) utilisation of a method 
for the re-projection of datasets and catalogues, which is very 
useful for the evaluation and registration of craters that are very 
close to poles; and (4) the extension of the previous framework 
for evaluation of CDAs (Salamuniccar et al., 2011b) with datasets 
and GT catalogue for the Moon, which can be used in parallel with 
those previously prepared for Mars. 

The new Crater Shape-based interpolation module improves 
our previous CDA as follows: (1) it decreases the number of false- 
detections for the required number of true detections, which 
leads to a smaller number of crater-candidates that need to be 
evaluated manually and at the same time to larger extensions of 
GT catalogues with previously uncatalogued craters; (2) it 
improves detection capabilities for very small craters, which leads 
to detection of considerably more craters from the same dataset; 
and (3) it provides more accurate automated measurements of 2D 
profiles and other crater properties such as depth-diameter ratio 
and more precise comparison of Lunar and Martian impact 
craters. In fact, the measure of relative difference in overall 
achieved performance between the CDAs from our previous work 
and this paper is 54.1% (wherein the only difference is in the new 
Crater Shape-based interpolation module - different from the 
linear one exclusively for distances smaller than one pixel), which 
is more than the difference in performance of 45.9% between the 
CDA from Stepinski and Urbach (2008) and our previous CDA 
(where the CDAs are completely different). A possible explanation 
- based on the presumption that CDA usually looks (in the first 
approximation) for circles/ellipses wherein all planetary bodies 
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Fig. 12. Depth/diameter in log/log scale for Martian (top-left) and Lunar craters (top-right), the average topographic-cross-profiles for 4 groups (centre), and more 
detailed insight (bottom). 


have craters with approximately those shapes - is that there is no 
reason why the same CDA should not work on different bodies. 
The experience of the other researchers is that this is not the case 
when the CDA is based on supervised learning, where the issue 
could also be the application of learning results from one image to 
another. This led to the investigation of advanced machine- 
learning techniques such as transfer-learning (Ding et al., 2011). 
The conclusion from this paper is that differences in shapes 
between Lunar and Martian craters (e.g. ~2 times larger/smaller 
average depth/diameter ratio) could lead to similar problems for 
other CDAs as well, because of e.g. the over-optimisation problem. 
From this perspective, instead of only claiming that the same CDA 
can be used for detection of craters from more than one planetary 
body, it is advisable to actually demonstrate this, as done in this 
paper for global DEMs of Mars and the Moon. 

The new Lunar crater catalogue from this paper is already in 
use in the current project of impact melt volume estimates of 
small- to medium-sized Lunar craters using LRO data (Mazarico 
et al., 2011). The Crater Shape-based interpolation CDA from this 


paper is already in use in current projects of cataloguing: 

(1) Phobos’ craters (Salamuniccar et al., 2011c, 201 Id); and 

(2) additional Lunar craters using DEMs reconstructed from 
optical images (Loncaric et al., 2011, Salamuniccar et al., 201 le). 
The goals and challenges for future work include integration of 
MA132843GT with the manually assembled catalogue, which 
currently contains over 288,993 craters expected to be completed 
soon (Robbins and Hynek, 2010), and the use of improved LOLA 
DEMs (1/256° and 1/512° currently, up to 1/1024° eventually) in 
combination with the newly released LRO global image mosaic, in 
order to improve the Lunar catalogue. 
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